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The bound state solutions of two fermions interacting by a scalar exchange are obtained in the 
framework of the explicitly covariant light-front dynamics. The stability with respect to cutoff 
of the J 7r =0 + and J^=f + states is studied. The solutions for J 7r =0 + are found to be stable for 

2 

coupling constants a = |— below the critical value a c (a 3.72 and unstable above it. The asymptotic 
behavior of the wave functions is found to follow a -pr+w l aw - The coefficient f3 and the critical 
coupling constant a c are calculated from an eigenvalue equation. The binding energies for the 
J' r =f + solutions diverge logarithmically with the cutoff for any value of the coupling constant. For 
a wide range of cutoff, the states with different angular momentum projections are weakly split. 

o 

£S) I PACS numbers: 11.80.Et,11.10.St,11.15.Tk 

I. INTRODUCTION 

One of the most difficult problems in field theory is the calculation of bound states due to the fact that they 

necessarily involve an infinite number of diagrams. A promising approach to deal with this problem is Light-Front 

Dynamics (LFD). In its standard version jlj, the state vector is defined on the surface t + z — 0. The bound states in 

ly-s i the Yukawa model (two fermions interacting by scalar exchange) were studied in references P, 0] using Tamm-Dancoff 

method (see, e.g., [0, pi). It was found in particular that, because the dominating kernel at large momenta tends 

04 i to a constant, the binding energy of the J = + state is cutoff dependent, what requires the renormalization of the 

r^~ Hamiltonian. 

The two-fermion wave functions were also considered in the explicitly covariant version of LFD (CLFD) [jj|. In 
this formalism, proposed in fit], the state vector is defined on the plane given by the invariant equation u-x — with 
oj 2 = 0. This approach keeps all along explicitly the dependence of the amplitudes on the light-front normal uj and 
i -£h ■ presents some advantages, in particular when calculating form factors H and when constructing non zero angular 
momentum states || . A first attempt to deal with CLFD wave functions was done in |l(j, [ll| where deuteron and pn 
j^-i. scattering J T = + state were calculated perturbatively and successfully applied to deuteron e.m. form factors (12] 
measured at TJNAF ||. 

The CLFD equations have now been solved exactly for a two-fermion system in the ladder approximation with 
different boson exchange couplings |L4[ |l5[| . The first results for the Yukawa model have been reported in |1(|. We 
investigated with special interest the stability of the bound state solutions relative to the cutoff, disregarding the 
self energy contribution and renormalization. We have found a critical phenomenon for the cutoff dependence of the 
binding energy. The J = + solutions were found to be stable - i.e. with finite limit when the cutoff tends to infinity 
- for coupling constant a below a critical value a c . On the contrary, for values exceeding a c , the system collapses. 
This fact manifests itself either as an infinite number of bound states with unbounded energies going to — oo for a 
finite value of the coupling constant or as a zero value of the coupling constant for a fixed value of the binding energy. 
The present paper is a detailed version of our work [[16| , includes new findings concerning the calculation of the 
critical coupling constant and the asymptotical behavior of the wave functions and the full treatment of the J = 1 + 
state. Our results are compared to those obtained in g] . 

In section II we remind the general properties of the two-fermion equations and wave functions in CLFD, already 
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presented in J6J. In sections [II and IV the system of equations for the J T = + and J"" = 1 + wave function 
components are derived. We show in section M that, after a linear transformation of these components and a change 
of variables, the CLFD equations are identical to the ones considered in |g|. In section VI we analyze analytically the 
asymptotical properties of the kernels and wave functions and their relation with the existence of a critical coupling 



constant. Numerical results are presented in section VII and section VIII contains the concluding remarks 



II. COVARIANT WAVE FUNCTION AND EQUATION 

We briefly describe here the main properties of the CLFD wave functions and equations. A detailed derivation can 
be found in §, [lj, [jj- 

In the covariant version of LFD the wave functions are the Fock components of the state vector defined on the 
light-front plane lj-x — 0. The standard LFD approach is recovered as a particular case with lj = (1, 0, 0,-1). 




FIG. 1: Graphical representation of the two-body wave function. 

The wave function < P a2 cr 1 of a two-fermion bound state is shown graphically in figure 0. It depends on four four- 
momenta 



^a 2 a t =$a 2 <?i(kl,k 2 ,P,OJT) (1) 

which are all on the corresponding mass shells (kf = k 2 = to 2 , p 2 = M 2 , (ujt) 2 = 0) and satisfy the conservation law: 

ki+k 2 = p + ljt. (2) 

In the standard approach the + and _L components are conserved, whereas the minus-component is not. These 
properties are reproduced by equation (g), since in this case the only nonzero component is oj_ ~ 2. Equation (g) is 
thus a covariant generalization of the usual conservation law. 

The general form of the wave function p) depends on the particular quantum number of the state and is obtained 
by constructing all possible spin structures. 

It is convenient to introduce variables (k,h), constructed from the initial four-momenta as follows: 



k = L^iV)^ =ki 
L- 1 {V)lj 



= [** 



kvP 



(3) 



where V = p + ljt, and L 1 (V) is the Lorentz boost into the reference system where V = 0. In these variables, the 
wave function (||) is represented as: 



*„ 



$« 



(k,n). 



(4) 



Under Lorentz transformations of the four-momenta ki,k2,p,0JT the variables k,n are only rotated J6|, so that 
parametrization (0), though being three-dimensional, is also explicitly covariant. In practice, instead of dealing 
with transformations (|3|), it is enough to write the wave functions and dynamical equations in the center of mass 
reference system, i.e. the one for which V — k\ + k 2 — and in which k\ = k, k 2 = —k, lj = n\£j\. 
The equation for the wave function in terms of variables (0) reads: 



M 2 - 4(fc 2 + to 2 )] <P a2ai (k,n) = —JY / K%% (k, k ', n, M 2 )<P a ,^ (k ', nf ' 

a', a' 



£k> 



(5) 



where K(k,k' ,n,M 2 ) is the interaction kernel. It depends on the scalar products of the vectors k',k,n and also on 
the scalar products k-a, k '-a and n-a. For the Yukawa model it will be precised in next section. 



In CLFD, the construction of states with definite angular momentum has some peculiarities which are explained 
in pi IT3. These peculiarities are related to the fact that the angular momentum operator J is not kinematical, 
but contains the interaction. However, we can overcome this difficulty by taking into account the so-called angular 
condition, derived from the transformation properties of the wave function under rotations of the light-front plane. 
Assuming that the state vector satisfies this condition, the problem results in finding the eigenfunctions of a purely 
kinematical operator M 
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-i k x 



d 


„ d 


- 


- i n x — 


dk 


an 



S\ + S 2 



(6) 



where si are the fermion spin operators. The operators (M 2 , M z ) have the same eigenvalues J(J + 1) and A than the 
full angular momentum operators (J 2 , J z ). 

As already mentioned, in any Lorentz transformation of the state vector, k and h undergo only rotations, with 
the same rotation operator than the one acting on spinor indices <7i,o"2. I n this respect, the eigenstates of M 2 ,M Z 
are constructed as in the non relativistic quantum mechanics. The only difference is that we have at our disposal 
two three-dimensional vectors (k,h) which enter in this construction on equal ground, instead of the only relative 
momentum k in the non relativistic case. 

Since the interaction kernel in (H) depends on the scalar products of all the three-vectors, including spin operator, 
M commutes with it. The solutions of the equation (||) are eigenfunctions of the operators M 2 , M z with eigenvalues 
J(J+1) and A. 

Although the operator M contains the derivatives both over d/dk and d/dn, its projection on n-axis does not 
involve any derivative with respect to variable n. Furthermore this latter enters in equation (|5|) as a vector parameter 
only and not as a dynamical variable. Therefore, there exists another operator which commutes with the kernel, 
namely: 

A 2 = (n-M) 2 . (7) 

Since A 2 is a scalar, it commutes also with M . Therefore, in addition to </, J z , the solutions are labeled by a: 

A 2 ip a (ic,n) = a 2 -ijj a (k,h). (8) 

The operator A 2 has J + 1 eigenvalues, a 2 — 0, 1, . . . , J 2 , and 2 J + 1 eigenfunctions which are split in two families 
of J and J + 1 states with opposite parities. For J = 0, there is only one value a — 0. For J 71 " = 1 + , there are two 
values a — and a = 1. The wave function with definite J is determined by Nj = \{2J + 1) x 2 x 2 = 2(2J + f ) 
spin components (the h factor comes from parity conservation). The system of equations for these 2(2J + 1) spin 
components is split in different subsystems each of them with definite value of a. For example, the wave function for 
J = 1 is determined by Nj = 6 components ]10| and the equation system is split in two subsystems corresponding to 
a = and a = 1 and containing respectively 2 and 4 equations (see section IV). 

The calculation technique of the CLFD is given by special graph rules which are a covariant generalization of 
the old fashioned perturbation theory. It was developed by Kadyshevsky |17| and adapted to CLFD in |g, 0. The 
equation for the wave function is shown graphically in figure and the corresponding analytical form (a) is obtained 
by applying the rules of the graph techniques to the diagram displayed in this figure. 

All along this paper we will consider the one boson exchange kernel with scalar coupling only. The interaction 
Lagrangian is C lnt = g ^ipcf)^' . The corresponding amplitude - represented graphically in figure has the form: 
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•■(k 2 )u^{k' 2 ) 



^(fciK^i: 



(9) 



*0) 



V(fci-M)) 



/i 2 -(fci-/ci) 2 + 2r'w(fci-fci) u 2 -(h - k[) 2 + 2TLo-(k 1 - k[) 
Like the wave function, the kernel is off energy shell for t,t' =/= 0. The two terms can be simplified into: 

^,2 



K(ki,k2,UT',k[,k' i 
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Am 2 (Q 2 + / j 1 2 ) 



[u CT2 (fc 2 )w^(fc2)] [u^kxju^ik'x)] 



(10) 



Q 2 is given in variables k, h by 

Q 2 = k 2 + k' 2 - 2kk' | 1 



-M 2 
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(gfc - £fc' 
2£fc£fc' 
kcos9 



cos 9 cos 9 1 — 2kk' sin 9 sin 9' cos ip' 
k' cos 9' 
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(11) 
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FIG. 2: Equation for the two-body wave function. 
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FIG. 3: One boson exchange kernel. 



where kcosO = n-k, k' cosO' — n-k' and tp' is the azimuthal angle between k and k' in the plane orthogonal to n. 



III. THE J = 0+ STATE 

The wave function of a two fermion system in the state J = + has the form jLl| : 

#<r 2 <7i (fel,&2 ,P,wt) = y/mu C r 2 (k 2 )4'U c u cri (ki), (12) 

with 

<t> = fiSi + fcS 2 , (13) 

/, are scalar functions depending on (k, 9) and S, are spin structures given by 



1 e k 2mui m 2 

m = 7= — 75, ^2 = F 9" 75 

2V2e k 2V2mksm6 \ w-p Ejg 



(14) 



in which u> — Uf/y^ and £& = \/k 2 + m? . U c = 7 2 7° is the charge conjugation matrix, u the usual spinor normalized 
as u a u a i = 2mS (7(7 i : 



1 



u a {k) = Vefc + m I g-k | u><r, 
(e-fc +m) 

and ui CT the two-component spinor normalized to u>J."U) CT ' = (W'- 



(15) 



Because of the spin structure u>, the relativistic wave function (1^) is determined by two scalar functions, in the 
agreement with the counting rule mentioned above. In section [V we will show that this state corresponds in the 
standard approach to the (1+, 2—) state discussed in B and described by the two components $ 1+ , $ 2_ . 



In the reference system where k\ + hz = 0, the wave function (12) can be represented as follows: 

with 

hit -\ l ( f i ^"[fcx"] t\ 



(16) 



The normalization condition has the form: 



El*- 



(2tt) 3 



d 3 k 
Tr{<j)(k2 + m)(j)(ki — m)} 

Sk 



Tr{ip* (fc, n)ip(k, n)} 



(2tt) 3 



(27T) 



JPk 



£k 



(17) 



where we denote k\ = fcf 7^ and <j> = 7o</>^7o- 

The spin structures Si in ( 14 ) are constructed in order to reproduce equation (jig ) without any additional coefficient 
and they are orthonormalized relative to the trace operation in (nT 



Tr{Si(k 2 + m)Sj(ki - to)} = S K 



where Si — 7o<S',- 7o: that is: 



S'i 



2^2, 
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w-p 



(18) 



(19) 



We insert expression (n3) for the wave function in equation (0) and multiply it on left by u(k2J and on right by 
u{k\). Using the relation J^ a w <T (fc)-u <T (fc) = k + m, we get: 



Af 2 -4(fc 2 +m 2 ) (fc 2 + m)0(fc 1 - to) 



(20) 



.T 



2tt 3 7 4to 2 (Q 2 + /! 2 ) 



(fc 2 + to)(/c 2 + m)4> (k[ — m)(ki — to) 



d 3 /fc' 



Substituting (/) in the form (|13|), multiplying equation (|o|) by Sj from (|l^) and using the orthogonality relation (|l§|), 
we obtain the system of two-dimensional integral equations for the components /;: 



M 2 -4(P + to 2 ) /i(fc,( 
M 2 -4(P + to 2 )] / 2 (fc,( 



to 2 /" d 3 k' 

£3 / [tf n (fc, 0; fc', 0')/i(fc', 0') + ^i 2 (fc, 0; k', e')h{k', 9')] ^- 

m 2 P d 3 k' 

—3 / [K 21 {k,e ] k , ,e , )h{k , ,6') + K 2 z{k,e-,k\e , )h{k , ,e')] — 

2tt 3 J ey 



The kernels Kij result from a first integration over the azimuthal angle ip' of more elementary quantities: 



Ka = 



1 



m'e k e k > J 



2t: 



Kij dip' 
f + m 2 2tt 



(21) 



(22) 



with 



k^ = ^■e fc e / fc rr[5 f -(fc 2 + m)(fc > 2 + m)S5(fc'i-m)(fci-m)J (23) 

However, for notation convenience, we keep in (El]) the three dimensional volume element though the kernels (K2) 
are ip' independent. The traces ( f23| ) are expressed through the scalar products of the available four- vectors. The 
analytical expressions for these scalar products and for Kij in the variables k, k! , 9, 9' , ip' are given in appendix El 



IV. THE J = 1+ STATE 

As mentioned in section II, the two-fermion wave function with J = 1 + is determined by six components associated 
to each of the six spin structures one can construct |Q : 

Q _ (fci - k 2 r „ i „ „ _ ^ (fci - ktYCj 

2to z m w-p Zmuj-p 

S B „ = ^l^ vpl k lv k 2p u~ n S 6fi = -. (24) 

The system of equations satisfying by these components can be split in two subsystems corresponding respectively 
to the eigenvalues a = 0, 1 of the operator A 2 (ph. The subsystem a = is - like in the J = case - determined by 
two components whereas the four remaining components are related to a = 1. The total number of components as 
well as the dimensions of the subsystems (2+4) coincides with the results obtained in the standard approach j|. 

A. The case a — 



In this section we consider the state with a — 0. We represent the wave function in a form similar to (12) 



^(Ai.fe.p.wr) = V^^(p,X)u^(k 2 )cj>^U c u^(k l ) , (25) 



e^{p, A) is the spin-1 polarization vector and functions (jr^ are represented as: 

- J°) o(°) , Jo) „(o) 



9rS%+g?>S%. (26) 

where g\ are scalar functions and the spin structures S\ will be expressed in terms of structures ( |24| ) . 
In the reference system where fci + fc 2 = this wave function can be represented in the form: 

<P x a2ai (k,h) = V^wl 2 ^\k,h)a y wt 1 . (27) 

One can easily check that the function ip°(k,n) - satisfying equation (pj) with a — - is proportional to n, i.e. it 
satisfies the condition ip° — n(n-ip°). Since we are dealing with a J" = 1 + pseudovector state, it has the following 
general decomposition: 

^n) = Jll^rt + ^Z'e^ V (28) 



in which the vector h is multiplied by the two only pseudoscalar structures one can construct a-k and a-h. 

The four-dimensional structures S\ 2 u, are Du 
from (Eq) written in the cm. -system. One finds 



The four-dimensional structures S\ J are built in terms of S^ defined in (|24|) in a way to obtain equation J28| ) 



l(0 ) V3M (0 ) v^M ^m 2 cos6>, 



s ^ ~ 27^ 53M S2 » ~ W^ \^k~ s ^ + 5e " ' ' (29) 



The normalization condition for the wave function reads: 



I E l^J-j^J^Tri^ik+rn^ih-m)}^ (30) 

A(72(Tl 



Tr{^(U)^°KU)}— = 7 ^ I (gf'f + ig^) 2 



3(2vr)3y ^ " e k (2tt)3 

where we have introduced the tensor 






n"" = I £ e "*(p> A ) e >' A ) - 5 (l? " 9 ^) 



(31) 



The structures S} are orthonormalized relative to the trace operation: 

WTr{S$(k + m)sfj{k - m)} = S %3 



(32) 



with Si — 7oSL 70 = S t . The same relation and the orthogonality condition (p3) hold for the structures SJ 
corresponding to a = 1 and constructed in the next section. 

We finally obtain for the components g\ a system of two equations having the same form than for J=0 - (|2l| ) and 
(E2n - with kernels Kij given by: 



Kij = 9 —e k e' k W v Tr sfKki + m)(k> 2 + m)S'^\k' 1 - m)(k - m) 



(33) 



Their explicit analytical expressions are given in appendix lAl 



B. The case a — 1 



The solution ip (k,n) corresponding to a — 1 is orthogonal to n, i.e. it satisfies n-ip 1 = 0. To fulfill this condition, 
it is convenient to introduce the vectors k± and <j± orthogonal to h: 



k — k cos Oh _ _, A A 

kj_ = - er_L = er - (n-cr)n 



and to write down ip 1 in the form: 



V3 r (i 



^(k, n) = ^ {sf^j. + .92 1} (2fe±(fc±-ax) - <7 ± ) + j^ fcj. (<?-n) + ^[fc x n]} 



The four-dimensional representation $„,— analogous to (129), is written in terms of 



,(D. 



A(l) - „(!) <?« + „(!) OW 4. „(!) <?« + „(!) <yU) 



The four spin structures 5* are orthonormalized according to (K52) and read: 



(34) 



(35) 



S il = Yl VijSjn, i = 1, ■ • • ,4; j = 1 



,...,6, 



(36) 



with Sjp defined in ( |24D and the coefficients r)ij given in appendix |B|. 

The normalization condition in terms of <pjj, and t/A 1 ) exactly coincides with ( p0| ) and in terms of the components 
<j^ is rewritten as: 



(20 



d 3 fc 



( 5 i 1) ) 2 + (3i 1) ) 2 + (^ ) ) 2 + (3i 1) ) 2 
The system of equations for the four scalar functions g\ has also the same form than (21) 



M 2 -4(/c 2 + m 2 ) 



£\k,n) = £j / ^K^(k,k',n)g^(k', 



d 3 k' 
) 

Eft' 



(37) 



(38) 



The kernels K\- are obtained from (33) by substituting S^ instead of S^°\ Their analytical expressions are given 
in appendix [A]. The details of calculation can be found in 113] . 

We are interested here in calculating the mass M of a J=l state. The physical solution satisfying the angular 
condition mentioned above is given by the superposition of the solutions (ip a , M a ) with definite a P, [15j: 



^=C O ^° + C 1 1 . 



The corresponding mass is 



2 J\ J? i „ 2 71^2 



M A = CqMJ + c\M{ 



V. RELATION WITH THE STANDARD APPROACH 

In reference B , the system of equations solved for the J z — helicity sector had the form: 
(M 2 -M 2 )$ x (R ± ,x) 

= ^ J [Vii{R^x\B!^x l )^ 1 {B!^x l ) + V 12 {R^x\B!^x l )^2{B! ± ,x l )]R! A _dBl A _dx l , 
(M 2 -M 2 )$ 2 (R ± ,x) (39) 

= ^ J [V 2l {R^x-R' ± ,x')^ 1 {R'^x') + V 22 {R^x-,R'^x')^ 2 {R' ± ,x')]R' ± dR' ± dx\ 

with Mq = -jrzrk ■ They correspond to equations (3.1a), (3.1b) from J2[ with the notations k = R±, q = R' ± , y = x'. 

The coupled wave functions ($i, $2) should be identified either to ($ 1+ , $ 2_ ) or to ($ 1 ~, $ 2+ ) sets. The kernels Vij, 
given by eqs. C1-C4 in Q, are for the case ($ 1+ , $ 2 ~): 

y. /" 27r %'<V 

y 7o x(l-.TK(l-.T')(Q 2 +/i 2 )' 
v u = i?±i?l(2xx' - a: - x') 

+ [R\x'{\ - x') + Rfx(l -x)-m 2 (x + x')(2 -x- x')] cos tp', 
v 12 = m [R±{x + 3a/ - 2xx' - 2a;' 2 ) - R' ± (x' + 3x - 2x'x - 2x 2 ) cosy'] , 
v 2 i = m \R'j_(x' + 3x — 2x x — 2x ) — R±(x + 'ix' — 2xx — 2x ) cosy'] , 
v 22 = R_lR' ± (2xx' - x - x 1 ) cos tp' (40) 

+ [R 2 j_x'{1 ~ x') + R'lx(l -x)- m 2 (x + x'){2 -x- x')} . 

The momentum transfer Q 2 reads: 

Q 2 = m 2 - (l - -) 2 + -R\ - 2R A _R' ± cosy' + ^R' 2 + (x' - x) ( "f* % - M 2 ) (41) 

x V x' I x x' \ x' ( 1 — X ) J 



for x < x' and with the replacements x «-> x' R±_ *-> R' ± for x > x' . 
The functions $, are normalized as: 



J (|*i (flx, x)\ 2 + \$ 2 (R±, x)\ 2 ) R± dR ± dx = 1. (42) 

We will show that the equations ( |2l|) with kernels (|Al[) corresponding to our J = state are identical to equations 
( p9| ) with kernels (pOJ) for ( < ! >1+ , $ 2_ ). They transform into each other by a linear combination of the wave functions 
and by a change of variable. The relation between the wave functions can be written in the form: 

f 3 = N{x)Y,*i<%(R±) ( 43 ) 



with a normalization factor ensuring the equivalence between the conditions (17) and (f42j) 

N(x) = 2 3/2 Trm- 1/2 ^x(l - x) 
and a unitary matrix c: 



1 (- R ± m \ (44) 

y/R 2 ± + to 2 V m ^-L / 



Variables (i?j_,a;) used in (J39|) are related to (fc,0) from (21) by 



1 / fccos^. 
R±=kwa.6, x = -[l (45) 

2 V £fc 



and their reversed relations: 



4x(l - ar) 



cos 9 



(l-2x)v^ 



v /i?i+m 2 (l-2x) 2 



(46) 



Variables (i?j_,a;) can be directly constructed from the initial four- momenta as follows. We introduce the four- 
vector R = k\ — xp with x = u}-ki/{ui-p) and represent its spatial part as R = R\\ + R±, where i?y is parallel to uj and 

R±-uJ = 0. Since by definition of R, Ruo — Rouiq — Rn-iv = 0, it follows that |i?o| = |-R|||, and hence R\ = —R 2 is 
invariant. 

For J = l,a = quantum numbers, corresponding to the state (1 _ ,2 + ) from J2|], the relation between the <r- 



components and $j is the same than (|43| 



(0) 



N{x)Y,c lJ (Ri_)^ l 



(47) 



with a c matrix 



k^R\ 



R±e k -k z m \ 
TT 7 \ —k z m —R±Ek J 



(48) 



Note that inJEj k z = —kcos9 with a minus sign due to h — —z). 

Inserting ( J43| ) into equations (pl|), we reproduce equations ( J39J ) with kernels Vy which are linear combinations of 



N(x') 

13 N{x)N 1 (x') ^ 
i'j' 



7 -^2ci:i{R±)Ki: j/ c j < j (R' ± ), N^x') = — x '(l-x') 



(49) 



Substituting here the kernels K^ from ( |22| ) for J = 0, we get the expressions for Vy given by (§{j|). A similar identity 
is obtained for the kernels Kij of J = l,o = state ( |A2|) and the kernels Vy- for the corresponding (1~,2 + ) state 
considered in B. Hence, both systems of equations are strictly equivalent. 



It is useful to explain the origin of the linear transformations (E3[) and (47). They result from using different 



representations of the spinors and different normalizations of the wave functions. Namely, instead of eq. ([15]) the 
following spinor is defined in [El : 



< t {h) = 



k 3 



mk + (fc + k 3 + 7o7 ± -fc- L )A^ 







where 



A ± = J 7 T 7 ± , 7*^70 ±73, k ± = (k 1 ,k 2 ). 



It can be transformed to the form: 



(50) 



<"(k) 



m + ko + o-k 1T3 



H ; n 



2(fc + fo) V (ko - rri) 0-3 + <?-fc 
In order to relate the two wave functions, one may express the spinors u LF in terms of u: 

u L /{k)= Y, /W'(fcK'(*0- 



(51) 



(52) 



With the explicit formulas (113) and (plj) we find 



1 



Paa'(k) = ^-U a >{k)u^ F {k) = „ 

2m y 2(fco + K3) 



ko +m + 
— wL 



d-ku-j, 

ko + m 



(53) 
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Note the opposite order of indices in both sides of ([53]). Below we give explicitly these matrix elements in the 
cm. -system k\ + fc 2 = 0: 

/? u (fci)=/3ii(fci)= vLZ™L =, /9i-i(fci)=-/?*4i(*i) = 



v /2fc+(fc +m)' 2 2 ~ 22 v /2fc+(fc +m)' 



(54) 

with fci + i/c 2 = k±e 1 ^ . Coefficients /3(fei) relate the spinors of particle 1 with momentum k\ = k, whereas coefficients 
/3(fe) relate the spinors of particle 2 with momentum fc 2 = —k). For the conjugated spinors u we get: 



w LF 



f(*i)= E /w*i)*m(*i). < F ( fc 2)= E / 3 ; 2 ^( fc 2)^(fc 2 ), (55) 



The matrix (3 is unitary: J^, /3j. .//9 (7 v" — So-' fit'vPa'o-" = &oo" ■ Substituting the spinors ( J55| ) into the wave function 
^ defined by means of the spinors u LF , we get the relation: 

* fflff2 = E #i*i (*i)/^oi (W„X (56) 

Note the opposite order of indices at \I> <Tl0 ' 2 and VV'o-' , i n accordance with the above definitions. The matrix elements 



of the wave function ip(k, ft) for J = 0, eq. (16), are the following 



^ii=Ci !=4s/2e~ 10 , Vi-i = ~^l-i = — 7=/i- (57) 

22 2 _ 2 V2 2 2 r 2 2 V2 



We substitute in (Kq) the coefficients (pi) and the wave function (57), following the paper @J, separate the phase 
factor: 

* <T1<J2 (i?- L , x) = E ^=$ CT1CT2 (i?,a;,m) 
V 2vr 

771 v 

and introduce the following linear combinations: 

<5> 1:t (R,x) = -=[*** (ii,s;-l)± $-5-3(^,^1)], 

v2 

$ 2± (i?,a;) = -=[$%-% (R,x;0)± $r^(R,x',Q)] (58) 

V2 



normalized according to (Q). In this way we reproduce for $ 1+ ,$ 2 the relation (43) with coefficients (|44|). For the 
other pair of components we get <I >1 ~ = <I> 2+ = 0. 

Similarly, taking the matrix elements of the wave function %p(°'(k,n) for J — 1, o = 0, equation (|28|), projected on 
the direction n (we remind that in the standard approach h\\ — z), we find: 

*m=*$ > l =^-( 9 ^R ± - 9 ^k 3 )e'^, 

2 2 2 2 y/2k 

2 2 2 2 ^/2fc 

In this way we reproduce for ($ 1_ , $ 2+ ) the relation J47| ) with coefficients (|4q), whereas we get <I >1+ = <I> 2 ~ = 0. 

VI. ASYMPTOTICAL PROPERTIES AND CRITICAL COUPLING CONSTANT 

The stability of solutions is determined by the asymptotical behavior of kernels Kij . This can be considered either 
for a fixed k(k') in the limit k'(k) — > oo or in the limit of both k, k' — > oo with a fixed ratio 7 = 4-. We illustrate in 
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what follows the analyzing method for the J = state. The relevant results for J — 1 are summarized at the end of 
the section. 



For a fixed value of k(k') we get: 



K^k, 0,^,0') 

K 12 (k,e,k',e') 
K 2 i(k,e,k' ) e') 

K 22 (k,9,k',9') 






if k 



'-12 

k 



C21 



k' fixed 



if k! — > oo, fc fixed 
if fc — » oo, fc' fixed 

if fc' — * oo, fc fixed 
if fc — » oo, kl fixed 

-#■ if fc' — ► oo, fc fixed 
c 22 if fc — > oo, k' fixed 

c 22 if A;' — * 00, k fixed 



(59) 



with coefficients en, C12, c 2 \, c 22 depending on k' , 9, 9' and c' 11; c' 12 , c 21 , c 22 depending on fc, 9, 9' . 
One has for instance 



airk' sin 6* sin 9' 



f'22 



m 2 (gfc' — fe' cos cos 0' + \ek' cos 9 — k' cos 0' | ) 



>0 



(60) 



and d 22 obtained form C22 by the replacement kl — ► fc,0 «-> 0'. Note that these behaviors are obtained once the 
integration over ip' in (E3) is performed. 

It follows from (pit) that the second iteration of K\\ converges at k' — ► 00: 



fCllGoifn 



rf 3 fc' 
£fc' 



1 1 1 fc'W 

Ffc^fc 7 fc' 



dfc' 



oc const. 



where Go ~ 1/k' 2 is the intermediate propagator. The integrals J K 2 iGqK\x^-^- , J K\\G§K\ 2 ^- and 
/ K 2 xGqK\ 2 ^-^- also converge whereas the second iteration of K 22 is logarithmically divergent: 



K 22 G K 22 - 



d 3 k' 



const ■ 



1 k n dk' 



/£-«,>. 



£fc' J A:' 2 k' 

The integral j Kx 2 GqK 22 ^-, J K 22 GqK 2 x^- and f Ki 2 GoK 2 i^-Z- are also logarithmically divergent, as a man- 
ifestation of the logarithmic divergence in the LFD box fermion diagram. 



In the domain where both k, k' tend to infinity with a fixed ratio j- — 7, it is useful to introduce the functions 



Ki„ 



7TCV 

in 



y/7An(e,#,i) if T < 1 
a\ -^^(5,^,1/7) if 7 >l 



where we have extracted for convenience the factor ^7. We find for An: 



A n (6,e',i) 



dip' 



27(1 — cos 9 cos 0') — (1 + 7 2 ) sin 9 sin 9' cos <// 



and for A 22 : 



A 9 



),e', 7 ) = - — 



1 
^/j J 27r (1 + 7 2 )(1 + |cos6»-cos6"| -cosflcosfl') - 27 sin sin 0' cos <// ' 

2w dip' (1 + 7 2 ) sin 6> sin 9' - 2 7 (1 - cos 9 cos 0') cos tp' 



1 



a/7 Jo 27r (1 + 7 2 )(1 + I cos0 — cos(9'| — cos0cos0') — 27 sin sin 0' cos y/ 



(61) 



(62) 



(63) 



with A n (9, 9', 7) ~ +^7 and A 22 (9, 9', 7) - -1/^7 for 7^0. 

Comparing the above formulas, we see that the dominating kernel is -K22 . It does not decrease in any direction of 
the (k, k') plane, whereas in the domain k — > 00, kl fixed (and vice versa) K\\ decreases. However, in the domain 
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^- = 7 fixed, none of the diagonal kernels decreases. The positive function ( J62| ) corresponds to an attractive kernel 
Kn whereas the unbounded function ( p3| ) correspond to a repulsive K22 (note that the relation (|61j) between iQj and 
Ajj contains the sign minus). 

The preceding results will be used to find the critical value of the coupling constant, above which the solutions are 
not stable as a function of the cut-off value k max . Since K22 is repulsive it cannot generate a collapse. The formalism 
is therefore developed in the one channel problem, e.g. the /1 component with the kernel K\\ satisfying the equation 



[M 2 - 4(fc 2 + rn 2 )] /(*, z) = ^ \ K(k, z; k' , z')f(k', z')-^. (64) 



d 3 k ' 
2n3J --,..,-,..,..„,..,..) — 

where z = cos 9 and the channel indices are hereafter omitted. 

Our further analysis leans on a method developed by Smirnov Jl^] and uses the fact that at k — > 00 the integral 
in r.h.-side of (p4|) is dominated by kl ex k. Indeed, when kl = fixed, k — > 00, the kernel K = K\\ in (p4T) decreases 
like 1/fc (see (|59|)), that results in the 1/fc 3 asymptotics for /. We will see below that when k! ex fc, the wave function 
decreases more slowly than 1/fc 3 : 

/(M~£^> 0</3<l. (65) 

Therefore, the integration domain kl ex k gives dominating contribution. Hence, in the integrand of eq. (B4T) one can 
replace the kernel by its asymptotics ( |6l| ) and substitute the asymptotical form of the wave function ( |65[ ) . 
Let us put k' — jk and consider the limit k — > 00 in which equation tends to 

to 2 r°° r +1 

-Af(k,z) = —J 7 d 7 y dz'K(k,znk,z')f( 7 k,z') 

where we have neglected the binding energy, supposing that it is finite. 

The integral over 7 can be split in two terms in order to use the limiting values of the kernel in its form ( |6l| ) 

a r 1 , , r r 1 , _„.„,..,,„ r , 1 



4/(fc,z) = -J dz'ij 7 d7x/7A(0,0',7)/(7* ; ,* , )+y i ^ — A(9,e',l/ 7 )f( 7 k,z') 
Both terms can be gathered by putting 7' = - in the second one 

4/(fc, z) = ^ y ^ dz' y d 7 A(5, 0', 7) \l 3/2 f{lK z') + 7-5/2/(^^1 (66) 

Inserting in (|6^) the wave function in the asymptotical form (|65| ) leads to the eigenvalue equation 

f +1 
h(z) = a / dz Hp(z, z) h(z') (67) 

with 

f 1 1 
H (z,z')= — — d<yA(z,z', 7 ) cosh(/31og 7 ) (68) 

io 2vr77 

The relation between the coupling constant a and the coefficient (3, determining the power law of the asymptotic 
wave function, can be found in practice by solving the eigenvalue equation ( p9| ) for a fixed value of j3 

,+1 
\ h(z)= dz'Hp(z,z')h{z') (69) 

and taking 

a(/3) = -L (70) 

A/3 



One should notice that the contribution 7 — v in the first term of integral (|66|) introduces corrections of higher order 
than i7rfw an d does not modify the asymptotics of the solution /. 
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The kernel A(z,z',j) in eq. ( p8[ ) is positive (see eq. (p2|)). The function cosh(/31og7) and, hence, the kernel 
Hp(z, z') have minimum at (3 = 0. This value of beta corresponds to maximal - critical - value of a (a c ). 
It is worth noticing that if A does not depend on z, z' , A(7, z, z') — ^(7), one gets [ |19| : 

i-/'*r^ Pi) 



(V, 



(i 



7T^/7 



Applying this condition to the non relativistic Schrodinger equation with potential 

a' a 1 
V ( r ) = — 2 = ~o 2' 

that corresponds to ^(7) = 1, one gets the critical coupling constant a c = ?, that is aj. = 7L- in agreement with 
p0[ . For the LFD Yukawa model, An is majorated by ^7. Inserting this value in (|7l|), we conclude that a c > n, in 
coincidence with estimation given in |lq| . 

Let us now examine the asymptotics of component f 2 . We consider equation (UW on a finite interval < k < k max 
and investigate the behavior of f 2 at k 3> m. Since kernel -^22, in contrast to -ftTn, tends to a constant at k — > 00 its 
contribution leads to the 1/fc 2 asymptotics for f 2 : 

f 2 (k,z)cc p . (72) 

It seems at first glance that with 1/k 2 asymptotics one gets a logarithmic divergence of the integral in the r.h.s. of 
equation (^lj) for f 2 : 

d 3 k , rk max 1 k n dk , 
K 2 2 h —~J o -JJ2 — lo g(^m- ) 



. However from the asymptotics of equation ( |2l| ) it follows that the asymptotical coefficient h 2 (k max , z) decreases as 

h 2 (k max ,z) 



l + b log(k r , 



thus ensuring the convergence of the r.h.s. integral in (21). 



For J = 1, a = state, the A 22 function behaves as ~ +-7= when 7 — > 0, what corresponds to a singular attractive 

-K22 kernel and generates a collapse. 

A similar situation occurs for J = 1, a = 1. Inspection of the four diagonal kernels shows that at 7 — > 0, function 
A 22 behaves also as A 22 ~ +-7^ namely: 

Ag= 1 ' a = 1 \9,0', 1 -,O)- l sin^sin 2 0' 



^4(1 + |cos0 - cos6*'| - cos 9 cos 9')' 
This singular attraction is responsible for the J — 1, a = 1 instability. 

VII. RESULTS 

Let us first present the results given by the single equation for /1 with kernel K\\ in the J = case. In all the 
calculations, the constituent masses were taken equal to m—1 and the mass of the exchanged scalar is /i=0.25. 

The numerical solution of equation ( |69| ) with the function An (7, z, z') given by ( |6l| ) is plotted in Figure |i| 

The critical coupling constant is obtained for (3 — for which the eigenvalue is A c = 0.269. It corresponds, according 
to (|7fj|), to a c = 3.72, in agreement with our numerical estimations a c > it fll6| . 

We have plotted in Figure || the mass square M 2 of the two fermion system as a function of the cutoff k max for 
two fixed values of the coupling constant below and above the critical value a c = 3.72. In our calculations the 
cutoff appears directly as the maximum value k max up to which the integrals in (|2l]) are performed. One can see 
two dramatically different behaviors depending on the value of the coupling constant a. For a — 3, the result is 
convergent. On the contrary, for a = 4, i.e. a > a c , the result is clearly divergent: M 2 decreases logarithmically as a 
function of k max and becomes even negative. This divergence is due only to the large k behavior of K\\. Though the 
negative values of M 2 are physically meaningless, they are formally allowed by the equations (^l|) and ([59|) . The first 
degree of M does not enter neither in the equation nor in the kernel, and M 2 crosses zero without any singularity. 
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4 r 




FIG. 4: Function a(/3) for LFD Yukawa Ku. 



The value of the critical a does not depend on the exchange mass fj,. For /i <§; to, e.g. fj, « 0.25, its existence is not 
relevant in describing physical states since any solution with positive M 2 , stable relative to cutoff, corresponds to 
a < a c . For fi ~ m one can reach the critical a for positive, though small values of M 2 . 

Let us consider now the two-channel problem. The kernel dominating in asymptotics is K22- In the case J = 
it is positive and corresponds to repulsion. Because of that, this kernel does not generate by itself any instability 
but cannot prevent from the collapse in the first channel (for enough large a), since due to coupling between the 
two channels the singular potential in channel 1 "pumps out" the wave function from channel 2. We would like to 
emphasize that the divergence in the J = case, when it happens, is not associated with the non decreasing behavior 
of the K22 kernel but with the existence of a critical value of the coupling constant separating two dynamical regimes. 
This property is due only to the attraction and large k behavior of K\\. 

In the coupled equations system ( |2l| ) the situation with the cutoff dependence is thus the same as for one channel. 
In Figure || is displayed the variation of M 2 for J = - or (1+, 2— ) and J = l,a = — or (1 — ,2+) states as a 
function of the cutoff k max . The value of the coupling constant is a = 1.184, the same that in Figure 2 of ||, below 
the critical value. Our numerical results are in agreement with those presented in this figure for a cutoff A < 100, but 
our calculation at larger k max leads to different conclusion for the J = state. One remarks a qualitatively different 
behavior of the two states. In what concerns J = 0, the numerical results become flatter when k max increases, with less 
than a 0.5% variation in M 2 when changing k max between k max =W and 300. The same kind of behavior is manifested 
when studying the cutoff dependence of the coupling constant for a fixed value of the binding energy. Figures and g 
show the a{k m ax) variation for B=0.05 and B=0.5. The J=0 state is very well fitted by a law a(k max ) = a^ + -gr 3 — 
with parameters a^ = 1.140, a = 0.0639 and v = 0.210 for B=0.05 and a^ = 2.648, a Q = 0.809 and v = 0.514*for 
B=0.50. We thus conclude to the stability of the state with J = 0, as expected from our analysis in section VI. 

We have examined the asymptotic behavior of the wave function and found that it accurately follows the power 
law (65) with a coefficient /3(a) given in Figure [|. For instance for a binding energy B=0.05 (a = 1.096) a direct 
measurement in the numerical solution plotted in Figure gives /3 = 0.820 ± 0.002 whereas the solution of equation 
([39]) for the corresponding a gives (3 — 0.819. The same kind of agreement was found for B=0.5 (a — 2.480): the 
asymptotic wave function - displayed in Figure JUl - gives /3 = 0.548 ± 0.002 and equation (p9) provides the value 
[3 = 0.547. This agreement shows in particular that the critical value of the coupling constant is the same for the 
one- and the two-channel problem. The influence of the second channel seems to have no any effect in the asymptotic 
behavior of f\. This channel behaves asymptotically as 1/fc 2 i.e. (3 — for any value of the binding energy, as 



indicated in sect. VI. One can see that the component fa changes the sign. 

It is worth noticing that - at least in the framework of this model - one could measure the coupling constant from 
the asymptotic behavior of the bound state wave function. We would like to point out however that the extraction of 
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M' 




FIG. 5: Cutoff dependence of the binding energy in the J = or (1+, 2—) state, in the one-channel problem (/i), for two fixed 
values of the coupling constant below and above the critical value. 



fc+m 



was 



coefficient j3 is numerically delicate when solving the equation with finite values of the cutoff k max . A way to overcome 

this difficulty consists in mapping the interval k € [0, oo] into a compact interval. The mapping k —* x 

used. 

Let us now consider the J = 1 case. For J = 1, o = or (1— ,2+) state, contrary to the J = 0, the value of 
M 2 (k max ) - displayed in Figure ^| - decreases faster than logarithmically and indicates a collapse. The asymptotics 
of the kernel K 22 ~ is the same than —K 22 ~ , it is negative and corresponds to attraction. The integral (|68|) for 



the kernel H(z, z') with the function A 22 given by (|63|) diverges. Therefore it results in a collapse for any value of 
the coupling constant, as pointed out in gj. The same result was found when solving the J — equations with the 

opposite sign of K 22 ~ ■ The case J = l,a = 1 state requires a four channel calculations. The results displayed in 



Figure 11 show a logarithmic divergence. 

One can see from Figures [7] and |ll| that the binding energies for the J = 1 states with different values of projection 
a are different but almost degenerate for a wide range of cutoff variation. For instance, with k max — 10 one has 
a a= o=1.17 a a= i=1.18 and with k max = 90, a a= o=l.lA a a= i = 1.16. These differences are less than 1% for a system 
with not so small binding energies (B=0.05) and we expect them to be negligible for weakly bound systems like 
deuteron. This quasi-degeneracy is much smaller in comparison to the results previously found with the scalar 
particles P, pll p3. In this latter case, a bound state with the same energy presents a splitting of rs 20% in a, what 
correspond to a energy difference AB « B. It is worth noticing that even for the sizeable cutoff k max = 90, the value 
of the coupling constant is still ~ 10% far from the converged one obtained by a mapping. 

The 2+4 components of the J = 1 + state are displayed in Figures [12 and O as functions of the momentum k . 
Solutions <7i and g 2 coming from the a = sector are comparable. They both depend on the angle 9 between the 
direction of the light-front plane and the momentum k: when k _L n, the first component vanishes and the second one 
reaches its maximal value as a function of 9: it dominates even over gi(k,9 — 0) and in the whole range of momentum. 
In the J = 1, a = state, only one solution is sizeable and dominates the three others for all values of 6. We would 
like to notice that the displayed components are related to the physical wave functions considered in Q by some 
specific linear combinations ensuring the correct transformation properties of the wave function, as well as an easy 
link with the non relativistic solutions. The peculiarities of these functions and their relations will be discussed in a 
forthcoming paper (la] . 
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FIG. 6: Cutoff dependence of M 2 for J = (l+,2- 
a = 1.184. 



and J = 1, a — (1 — ,2+) states, in full (two-channel) problem, for 



VIII. CONCLUSION 



The Light-Front solutions of the two fermion system interacting via a scalar exchange have been obtained. We have 
found that the J = - or (1+, 2—) - state is stable (i.e. convergent relative to the cutoff k max — > oo) for coupling- 
constant below some critical value, in a way similar to what is known in non relativistic quantum mechanics for the 
—a'/r 2 potential. In this point, our conclusion differs from the one settled in B, where it was stated that the integrals 
in eqs. (p9) diverge logarithmically with cutoff. Above the critical value the system collapses. This fact manifests as 
an unbounded value of M 2 when the cutoff tends to infinity. 

We have shown analytically that the asymptotic behavior of the wave function has the form , 2+/) and that the 
relation between coefficient /3 and the coupling constant can be obtained as a solution of an eigenvalue equation 
suggested by |0|. This relation provides in particular the critical value of the coupling constant, which corresponds 
to (3 — 0. These results are in agreement with the numerical solutions of the LFD equations. 

In the J = 1, a = - or (1—, 2+) - state the system is found to be always unstable, as it was pointed out in [||. 
The instability is related to the dominating K22 kernel which is attractive. The origin of the collapse is thus different 
from J = state, for which the K^ 2 ~ kernel is repulsive and the instability is due to the asymptotic behavior of 
attractive K\\ and depends on the value of a relative to a c . 

The solutions for the four channel problem J = l,a = 1 state have been obtained. They show also a logarithmic 
divergence of M 2 (k max ) and are thus unstable without regularization. 

These results should be taken into account when carrying out the renormalization procedure. The explicitly 
covariant version of Light-Front Dynamics (CLFD) seems very promising for handling this problem p3| , like it has 
proved to be fruitful in the Yukawa model. 
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FIG. 7: Cutoff dependence of the coupling constant, for FIG. 8: Cutoff dependence of the coupling constant, for 
J = and J — l,o = states, in full (two-channel) J = and J = l,a — states, in full (two-channel) 
problem, for B = 0.05. problem, for B = 0.5. 
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FIG. 10: Asymptotic behavior of the J=0 wave function 
components /, for B=0.5, a=2.48, /i=0.25. The slope 
coefficient are fii = 0.55 and P2 ~ 0. 



Joseph Fourier, in Grenoble, where this work was performed. This work was partially supported by the grant No. 
99-02-17263 of the Russian Fund for Basic Researches. The numerical calculations were performed at CGCV (CEA 
Grenoble) and IDRIS (CNRS). We thank the staff members of these organizations for their constant support. 



18 



1.185 



a 




FIG. 11: Cutoff dependence of the coupling constant for J 
logarithmic divergence. 



1, a 



1 state with B — 0.05 until k„ 



=300. It shows a 



APPENDIX A: EXPLICIT EXPRESSIONS FOR LFD KERNELS 



The kernels are obtained through the traces (|23|). To calculate them, we express the scalar products between the 
different four-vectors in terms of variables k, k' , 9, 9', ip' . The following results are useful: 



u z = 

,2 



n ■ i — n 9 — n. i — rv 9 — Til , 

LO-kl — x u>-p, 
uj-k 2 = (1 — x)to-p, 
uj-ki = X u>-p, 
oj-k' 2 = (1 — x')u)-p, 
kyki = 2e\ - m 2 , 

12 ~~ k' ' ^^ 1 

k v p = 2^.(1 - x) + M 2 x/2, 

k 2 -p = 2e 2 k x + M 2 (l-x)/2 

k\-p = 24' (1 - x) + MV/2, 

^•p = 24,a;' + Af 2 (l-a;')/2, 

fci-fei = -fcfc'sinflsinfl'cos^' + 2e^x + 2e^x' - 2e^xa;' - 2e\,xx' , 

k 2 -k! 2 = -fcfc'sinflsinfl'cos^' + 2e^x + 2e^,x' - 2e\xx' - 2e\,xx' ', 

k y k 2 = kk' sin 9 sin 9' cos <p' + 24(1 - ir)(l - x') + 2e 2 k ,xx', 

k 2 -k[ = kk' sin 9 sin 6>' cos <p' + 24,(1 - z)(l - x') + 2e 2 k xx\ 

where 2x = ( 1 — fccos6> J and similarly for x' . 

The analytical expressions for the Kij kernels (g3|) in the J = case are: 

K11 = -air [2fc 2 fc' 2 + 3fc 2 m 2 + 3fc' 2 m 2 + 4m 4 - 2kk' e k e w cos 6> cos 9' - kk'(e 2 k + £%,) sin6>sin6»' cos<//] 
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FIG. 12: Components g\ and <?2 of the J — I, a = FIG. 13: The four components gi of the J = 1 



= 1 



(1-, 2+) state, for B = 0.05 



(l+,2-) state, for B = 0.05. 



K12 = — anm(k 2 — k' 2 ) (k' sin9' + fc sin # cosy') 

«2i = -airm(k' 2 - k 2 ) {ksm.0 + k' sin 9' cos tp') (Al) 

-cutt [(2fc 2 fc' 2 + 3fc 2 m 2 + 3fc' 2 m 2 + 4m 4 - 2kk'e k e k - cos6»cos6»') cosy' - kk'(e 2 k + e\, ) sin sin 6'] 



«22 



where we have used a = §-. 

The analytical expressions for the Kij kernels ( p3|) in the J = 1 + , a = case are: 

Kll = - a7 r { [-2fcfc'e fe e fc / + (2fc 2 fc' 2 + 3m 2 k 2 + 3m 2 k' 2 + 4m 4 ) cos 6» cos 6»'] + e k e v (k 2 + k' 2 + Am 2 ) sin 6» sin 6»' cos y' } 

K12 = -cnrm [-£fc'(3fc 2 + fc' 2 + 4m 2 )cos6»sin6>' + £ fc (fc 2 + 3fc' 2 + 4m 2 ) sin cos 0' cos y'] 

k 21 = -ayrm [-e fe (fc 2 + 3fc' 2 + 4r7j 2 )sin0cos0' + e fe /(3fc 2 + k' 2 + 4m 2 ) cos 9 sin 8' cos tp'] (A2) 

K22 = -a7r{e fc £ fc /(fc 2 + fc' 2 + 4m 2 ) sin sin 0' + [-2kk'e k e k > + (2k 2 k' 2 + 3m 2 k 2 +3m 2 k' 2 + Am 4 ) cos 9 cos 9'] cosy'} 

The analytical expressions for the Kij, determining by eq. ( |22| ) the kernels K^ in the J — 1 + , a = 1 case are: 

2«u 



2«i 



= -m[e fe sin 2 6'(£ 2 . + 3e 2 .,) + e k > sin 2 0'(3ef +e 2 /)] - (cos 2 9 + cos 2 9' ){m 2 (e 2 k + e 2 k ,) + 2e 2 k e 2 k ,] + 4e k s k *kk' cosdcosd' 

— (e k + e k ') 2 sin $ sin 6' cos y' {(e^ — m)(ey — m) [cos cos 6*' + sin sin 0' cos 99'] — fc/c'} 

= me k (e\ + 3e 2 k ,) sin 2 (9 - me k -(ie\ + e\,) sin 2 0' + (cos 2 - cos 2 9')\m 2 e\ + m 2 e\, + 2e\e 2 k ,\ 

— sin0sin0' {kk'(e k — e k ') 2 + (e k + e k ') 2 (e k ' — m)(e k ' — m) cos0cos0} cosy' 



sm f(£fc — m 



2 k 



) {-(efc' + m)(£ fe - e fc ') 2 + <£fc' - m )(^fe + £fc') 2 cos 2 9'} cos 2 y 



13 



V2c 



2m 



= —sin 9' [2e k e k 'kk' cos 9 + (ey — m)[m(e\, + e k ) — 2els k <] cos0'] 

+ sin0 {2e k e k ikk' cos 9' + (e k — m)cos8 [e k '(el + e\, — 2me k ) — (e k > — m)(e k + e k i) 2 cos 2 6'] } cosy 
— (e k + e k ') 2 (e k — m)(e k > — m) sin 2 9 sin 9' cos 9' cos 2 ip' 



rnr 



= V2(k 2 - k' 2 ) {(e fc sin 2 6 + mcos 2 9)k' sin 9' + mk sin 8 cos tp' - k'(e k - m) sin 2 9 sin 9' cos 2 tp' } 

= -me k sin 2 6»(3e 2 , + e 2 ) + me k > sin 2 6»'(e^ + 3e|) - (cos 2 6» - cos 2 9') (m 2 (e\, + e 2 k ) + 2e^e|) 

— sin 9 sin 0' {fc/c'(efe — e k >) 2 + (e k + e k i) 2 (e k i — m)(e k — m) cos 8 cos 9'} cosy' 

+ sin 2 9'(e k ' - m) {-(e k > - e k ) 2 (e k + m) + cos 2 6»(e fc ' + e k ) 2 (e k - m)} cos 2 y' 



20 



2k 

— m(ek + ek>) 3 ~ (£k -m)(m(sl+e 2 k ,) -2sksl>) cos 2 9 - (ek< -m)(m(el+ el,) -2elek<) cos 2 9' 



an 

^2 



— Akk'skEk 1 cos 9 cos 9' — (ef. + £*/) {{e k — m)(e k ' — to) — kk' cos 9 cos 9'} sin 6* sin 9' coscp' 

- {sin 2 0' (sin 2 9e k e k ,(s 2 k + e 2 ., + 2m 2 ) + me fe ,(cos 2 9 + l){3e 2 k + e 2 k ,)) + (cos 2 9' + 1) 

x (me k sin 2 9(e 2 k + 3e 2 ,,) + (cos 2 9 + l)(m 2 (e 2 k + e 2 ,) + 2e 2 .£ 2 ,)) - %kk'e k e k > cos cos 0'} cos 2 y/ 



2k 23 



{2£fe£fe'fc'fccos6' + (efc/ — m)cos9'(ej,,m + £ 2 (m — 2£fe/))} 



Sin a — sm I 



\>2air 

x { (£fc — m) cos 6 ((gfc + £/t' ) 2 (£fc' — to) cos 2 0' — £fc' (e 2 + £ 2 / — 2mE k )) — 2kk'e k s k > cos 9' } cos <// 

— sin0' {(e k > — m)cos9' ((e k — e k >) 2 (e k + to) — (e k + e k ') 2 (e k — to) cos 2 0) + Akk'e k e k > cos0} cos 2 <// 

= \2(e 2 k , — £ 2 ) { mk sin 9 cos <// — fc' sin 9' \(e k sin 2 9 + m cos 2 0) — [(e/. +to) — (e k — to) cos 2 0] cos 2 </?'] } 

2k 3 i 



sin 9 [(ek — to) (toe 2 + me 2 , — 2e k e\, ) + 2sk£k' kk' cos ( 



\/2air 

+ sin 0' { (e k - to) cos 0' [(me 2 , + me 2 , - 2e 2 ,e fc / ) cos 2 9 + e k (e 2 , + e\, - 2me k , ) sin 2 9] + 2e k e k , kk' cos 9' } cos ip' 
— (e k + £fc' ) 2 ( £ fe — m ) ( £ fc' — to) sin 9 cos sin 2 9' cos 2 <// 



= sin0 [2£ fc £fc'fcfc'cos#' + (e — m)[m(e 2 k + e\,) — 2e k e 2 k \ cos 



\f2a-K 

+ sin 9' cos ip' {(e k ' — to) cos 9' [cos 2 9(me\ + me\, — 2e 2 e 2 ,,) + e k sin 2 9\e\ + e\, — 2ms k '] + 2e k e k 'kk' cos 9} 

+ sin0cos 2 tp \{ek — to) cos9 [ek'(2ekm — £ 2 , — e 2 ,) sin 2 9' + (2e k s\, — me\ — ms\,)(\ + cos 2 9')] — 4e k s k >kk' cos 9'} 



OTT 



= J£fc(£fc' - to) (e 2 +£ fe '(£ fe ' - 2m)) cos 2 0' - e k >m(3el + £ 2 ,) cos 2 - e k e k '(e 2 k + £ 2 , + 2m 2 ) sin 2 
((£fe + £fc') 2 (£fe — ni)(e k i — to) cos cos 0' — 2kk'e k e k ,) (cos cos 0' + sin sin 0' cos <//)} cos 93' 



= (e 2 - £ 2 / ) (e k — m) cos 9 sin 0fc' sin 9' sin 2 (/?' 

— il = V2(e 2 k , - £ 2 ) {mfc'sin0'cos</ + fcsin0 [(e fe / sin 2 0' + mcos 2 0') - (e fc / - m) sin 2 0' sin 2 <p'~\ } 

-^ = V^e 2 . -£?0 {mfc'sin^'cos^' + fcsinfl [[(m + £ fc /) + (to - s k > ) cos 2 6»'l cos 2 w' - (s k > sm 2 9' + mcos 2 9')]} 
a-K 

— L = — kiei — £?/)(£fe' — m) cos 6*' sin sin 0' sin 2 u/ 
a-K 

= {( £ fe + £ 2 ')(kk' sin 9 sin 9' cos ip' — m 2 ) + 2ek£k'(kk' cos 9 cos 9' — £&£&')} cos^p' 

APPENDIX B: COEFFICIENTS rfo DETERMINING THE ORTHONORMALIZED SPIN STRUCTURES 

The construction method of the four orthonormalized spin structures S\ ( p6|) determining the wave function cpy, 
(pq), is explained in p4[ . These structures are expressed in terms of six structures Sj^, in the form (p6|). The non 
zero coefficients rjij (i = 1, . . . , 4; j = 1, . . . , 6) are given below (z — cos 9): 

Vim 2 V3m V3(e k - ra)Ue 2 k + M 2 )z 

r ' 11 = ~A — 1 i T' ^ 12 = "1 ' 7 ? 13 = t« 2; ' 

4£fc(£fc + m) 4£ fe 16£^,fc 

= V3(4£ 2 +Af 2 ) ^ 

8£fem 

\/3TO 2 [£ fc (l-z 2 )+TO(l + z 2 )] V3; 

%l = : — TTTI ^ ) ^22 = 



//•/ 



4£ fe fc2(l-z2) ' i" 4s k 



V3(4£ 2 + M 2 )[ £fc (l - z 2 ) + m(l + z 2 )]z ^ v^r, 

16£ 2 fc(l- 
V3(4£ 2 +Af 2 )(l + z 2 ) 



7723 16£ 2 fc(l-z 2 ) ' ^ fc(l-z 2 ) : 



8£fcTO(l — z 2 ) 
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V3m 2 z V3(sk - m)(4el + M 2 )z 2 

mi 

7? 45 = - - y ^Z= = (BD 



V3m 2 z 




V3{e k - 

V3(4el + M 2 )z 
4£ fc ™^2(l - z 2 ) 


-m){Asl 
V2(l- 


+ M 2 )z 2 


2e k (s k + m)y/2(l 
\/3m 


-z 2 ) 

me = 


-z 2 ) 


k^2{\-z 2 )' 
a/3to 2 




2e k k^/2(l - z 2 ) 
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